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f_\ ' Abstract. We study the nucleon-nucleon (NN) problem as a six-quark system in 

^^ . a nonrelativistic chiral constituent quark model where the Hamiltonian contains a 

linear confinement and a pseudoscalar meson (Goldstone boson) exchange interaction 

between the quarks. This interaction has a long range Yukawa- type part, depending 

on the mass of the exchanged meson and a short range part, mainly responsible for the 

good description of the baryon spectra. We calculate the NN potential in the adiabatic 

^ ' approximation as a function of Z , the separation distance between the centres of the 

y^ . two three-quark clusters. The orbital part of the six-quark states is constructed either 

from the usual cluster model states or from molecular orbital single particle states. 

The latter are more realistic, having proper axially and reflectionally symmetries. In 

both cases the potential presents an important hard core at short distances, explained 

Q^ . through the dominance of the [51]i?5 configuration. However in the molecular orbital 

0^ ' basis the core is less repulsive, as a consequence of the fact that this basis gives a better 

upper bound for the energy of the six-quark system. We calculate the potential for 

the ^5'i and ^Sq channels with two different parametrizations. We find a small (few 

MeV) attractive pocket for one of these parametrizations. A middle range attraction 

^ ' is simulated by the addition of a cr-meson exchange interaction between quarks, of 

P5 . a form similar to that of the pseudoscalar meson exchange. The present study is an 

intermediate, useful step towards dynamical calculations based on the resonating group 

method. 



INTRODUCTION 

There have been many attempts to study the nucleon-nucleon interaction starting 
from a system of six interacting quarks described by a constituent quark model. 
These models explain the short range repulsion as due to the colour magnetic part 
of the one gluon exchange (OGE) interaction between quarks and due to quark 
interchanges between two 3g clusters [1,2]. To the OGE interaction it was necessary 
to add a scalar and a pseudoscalar meson exchange interaction between quarks of 
different 3g clusters in order to explain the intermediate- and long-range attraction 
between two nucleons [3-5]. 

^^ Talk presented by D. BARTZ at the International Workshop on Hadron Physics : "Effective 
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In [6] we have calculated the nucleon-nucleon (NN) interaction potential at zero- 
separation distance between two three-quark clusters in the frame of a constituent 
quark model [7-9] where the quarks interact via pseudoscalar meson i.e. Goldstone 
boson exchange (GBE) instead of OGE. An important motivation in using the 
GBE model is that it describes well the baryon spectra. In particular, it correctly 
reproduces the order of positive and negative parity states both for nonstrange [8] 
and strange [9] baryons where the OGE model has failed. 

The underlying symmetry of the GBE model is related to the flavour-spin 
SUf{S) X SUs{2) group. Combining it with the S^ symmetry, a thorough anal- 
ysis performed for the L = 1 baryons [10] has shown that the chiral quark picture 
leads to more satisfactory fits to the observed baryon spectrum than the OGE 
models. 

The one-pion exchange potential between quarks appears naturally as an itera- 
tion of the instanton induced interaction in the t-channel [11]. The meson exchange 
picture is also supported by explicit QCD latice calculations [12]. 

Another motivation in using the GBE model is that the exchange interaction 
contains the basic ingredients required by the A^A^ problem. Its long-range part, 
required to provide the long-range A^A^ interaction, is a Yukawa-type potential 
depending on the mass of the exchange meson. Its short-range part, of opposite 
sign to the long-range one, is mainly responsible for the good description of the 
baryon spectra [7-9] and also induces a short-range repulsion in the A^A^ system, 
both in the ^5*1 and ^5*0 channels [13]. This study is an extention of [6] and we 
calculate here the interaction potential between two 3g clusters as a function of 
Z, the separation distance between the centres of the clusters. This separation 
distance is a good approximation of the Jacobi relative coordinate between the two 
clusters. Under this assumption, here we calculate the interaction potential in the 
adiabatic (Born-Oppenheimer) approximation, as explained below. 

A common issue in solving the A^A^ problem is the construction of adequate six- 
quark basis states. The usual choice is a cluster model basis [1,2,14]. In calculating 
the potential at zero-separation distance, in Ref. [6] we used molecular-type orbitals 
[15] and compared the results with those based on cluster model single-particle 
states. The molecular orbitals have the proper axially and refiectionally symmetries 
and can be constructed from appropriate combinations of two-centre Gaussians. By 
using molecular orbitals, in Ref. [6] we found that the height of the repulsion reduces 
by about 22 % and 25 % in the ^5*1 and ^5*0 channels respectively with respect to 
cluster model results. It is therefore useful to analyse the role of molecular orbitals 
at distances Z ^ 0. By construction, aX Z ^ oo the molecular orbital states are 
simple parity conserving linear combinations of cluster model states. Their role is 
expected to be important at short range at least. They also have the advantage 
of forming an orthogonal and complete basis while the cluster model (two-centre) 
states are not orthogonal and are overcomplete. For this reason we found that in 
practice they are more convenient to be used than the cluster model basis, where 
one must carefully [14] consider the limit Z ^ 0. Here too, for the purpose of 
comparison we perform calculations both in the cluster and the molecular model. 



In the following section we recall the procedure of constructing molecular orbital 
single-particle states starting from the two-centre Gaussians used in the cluster 
model calculations. Then the GBE Hamiltonian is presented. The subsequent 
section is devoted to results obtained for the A^A^ potential. Next we introduce 
a middle range attraction through a scalar meson exchange interaction between 
quarks parametrized consistently with the pseudoscalar meson exchange. The last 
section is devoted to a summary and conclusions. 

SINGLE-PARTICLE ORBITALS 

In the cluster model one can define states which in the limit of large intercluster 
separation Z are right R and left L states 



R = -ibxr and 





In the simplest cluster model basis these are ground state harmonic oscillator wave 
functions centered at Zjl and —Z/2 respectively. They contain a parameter (3 
which is fixed variationally to minimize the nucleon mass described as a 3g cluster 
within a given Hamiltonian. The states (1) are normalized but are not orthogonal 
at finite Z. They have good parity about their centers but not about their common 
center r = 0. 

From R and L one constructs six-quark states of given orbital symmetry [/]o- 
The totally antisymmetric six-quark states also contain a flavour-spin part of sym- 
metry [/]f5 and a colour part of symmetry [222] c. In the cluster model the most 
important basis states [13] for the Hamiltonian described in the following section 
are 

R'L^ [6]o [33]^^) (2) 

R'L' [42]^ [33]^^) (3) 

R'L^ [42]^ [51]^^) (4) 

R'L^ [42]o [411]^^) (5) 

Harvey [14] has shown that with a proper normalization the symmetry [6]o contains 
only s^ and [42] o only s'^p^ configurations in the limit Z ^ 0. 

According to Ref. [15] let us consider also molecular orbital single-particle states. 
Most generally these are eigenstates of a Hamiltonian Hq having axial and refiec- 
tional symmetries characteristic to the A^A^ problem. These eigenstates have there- 
fore good parity and good angular momentum projection. As in the cluster model 
basis where one uses the two lowest states R and L, in the molecular orbital ba- 
sis we also consider the two lowest states, a, of positive parity and tt, of negative 
parity. From these we can construct pseudo-right r and pseudo-left / states as 



2"^/2 (a ± tt) for all Z, 



(6) 



where 



< r\r >=< l\l >= 1, < r\l >= 0. 



(7) 

In principle one can obtain molecular orbital single particle states from mean field 
calculations (see for example [16]). Here we approximate them by good parity, 
orthonormal states constructed from the cluster model states fl) as 



a 
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[2{1± < R\L >)]-'/ '{R ± L), 



Such molecular orbitals are a very good approximation to the exact eigenstates of 
a "two-centre" oscillator frequently used in nuclear physics or occasionally [17] in 
the calculation of the A^A^ potential. They provide a convenient basis for the first 
step calculations based on the adiabatic approximation as described below. 
Introduced in (6) they give 

R+L R~L 



r 
I 



± 



l+<i?|L>)i/2 (1- < i?|L >)i/2 



(9) 



At Z ^ one has a ^ s and tt ^ p (with m = 0, ±1) 
oscillator states. Thus in the limit Z ^ one has 



r 



2'/\s±p), 



where s and p are harmonic 



(10) 



and at Z ^ oo one recovers the cluster model basis because r -^ R and i ^ L. 

From (r, /) as well as from {a, it) orbitals one can construct six-quark states 
of required permutation symmetry, as shown in Ref. [15]. In the limit Z — ^ 
six-quark states obtained from molecular orbitals contain configurations of type 
gUpG-n -^j^j-^ n = 0,1, ...,6. For example the [6]o state contains s^, s^p"^, s^p'^ and 
p^ configurations and the [42] o state associated to the ^'-channel contains s'^p^ 
and s^p'^ configurations. This is in contrast to the cluster model basis where [6]o 
contains only s^ and [42]o only s'^p^ configurations, as mentioned above. 

Besides being poorer in s'^p^^'^ configurations, the number of basis states is 
smaller in the cluster model although we deal with the same [/]o and [/]f5 sym- 
metries and the same harmonic oscillator states s and p in both cases. This is due 
to the existence of three-quark clusters only in the cluster model states, while the 
molecular basis also allows configurations with five quarks to the left and one to 
the right, or vice versa, or four quarks to the left and two to the right or vice versa 
(see Eqs (11)-(19) of [18]). At large separations these states act as "hidden colour" 
states but at short- and medium-range separation distances they are expected to 
bring a significant contribution, as we shall see below. The "hidden colour" are 
states where a 3g cluster in an s^ configuration is a colour octet, in contrast to the 
nucleon which is a colour singlet. Their role is important at short separations but 
it vanishes at large ones (see e.g. [14]). 



HAMILTONIAN 

The GBE Hamiltonian considered in this study has the form [8,9] : 

H = Y.^^ + Y.§-.- ^^ + E Konf(r.,) + E Vxin,) , (11) 

with the hnear confining interaction : 

Konflr,,) = -^A,^ ■ A^^ {Vo + Cn, ), (12) 

and the spin-spin component of the GBE interaction in its SUf{3) form : 

{f=i 

+ E VK{n,)XfXf + V,{n,)X^X'^ + V,in,)X^X^] a, ■ a„ (13) 



F=4 



with A*^ = \/2/3 1, where 1 is the 3x3 unit matrix. The interaction (13) contains 
'-f = 7r,K,ri and rj' meson-exchange terms and the form of V^ {rij) is given as the 
sum of two distinct contributions : a Yukawa-type potential containing the mass 
of the exchanged meson and a short-range contribution of opposite sign, the role 
of which is crucial in baryon spectroscopy. 

In the parametrization of Ref. [8] the exchange potential due to a meson 7 has 
the form 

"■('') = ilai^f^C - ''""'? 9^ - ;^"^-p(-"V - r.?)}. (14) 

The shifted Gaussian of Eq. (14) results from a pure phenomeno logical fit (see 
below) of the baryon spectrum with 

ro = 0.43/m, a = 2.91 /m~\ (15) 

For the parametrization of Ref. [9], the potential has the form 

where 

A^ = Ao + K /i^, Ao = 5.82 fm-\ k = 1.34. (17) 

In the following, we shall call the form (14) Model I and the form (16) Model II. 
For a system of u and d quarks only, as it is the case here, the /^-exchange does not 
contribute. The apriori determined parameters of the GBE model are the masses 



TABLE 1. Parameters of the Hamiltonian (11-18) 





Model 


Vo (McV) 


C (fm-^) 


g'i/^^ 


.9o/47r 


Ref. 


I 





0.474 


0.67 


1.206 


8 


II 


-112 


0.77 


1.24 


2.765 


9 



mu,d = 340 MeV, /i^ = 139 MeV, /i^ = 547 MeV, n^> = 958 MeV. (18) 

The other parameters are given in Table I. 

It is useful to comment on Eqs. (14) and (16). The coupling of pseudoscalar 
mesons to quarks (or nucleons) gives rise to a two-body interaction potential which 
contains a Yukawa-type term and a contact term of opposite sign (see e.g. [19]). The 
second term of (14) or (16) stems from the contact term, regularized with param- 
eters fixed phenomenologically. Certainly more fundamental studies are required 
to understand this second term and attempts are being made in this direction. 
The instanton liquid model of the vacuum (for a review see [20]) implies point-like 
quark-quark interactions. To obtain a realistic description of the hyperfine interac- 
tion this interaction has to be iterated in the t-channel [11]. The t-channel iteration 
admits a meson exchange interpretation [21]. 



RESULTS 



We diagonalize the Hamiltonian (11)-(18) in the six-quark cluster model basis 
and in the six-quark molecular orbital basis for values of the separation distance Z 
up to 2.5 fm. Using in each case the lowest eigenvalue, denoted by {H)z we define 
the NN interaction potential in the adiabatic (Born-Oppenheimer) approximation 



as 



V^r^{Z) = {H)z-2mN-K, 



rel 



(19) 



Here rriN is the nucleon mass obtained as a variational s^ solution for a 3g sys- 



;ii)- 



oc 



tem described by the Hamiltonian 

exp [— (p^ + A^) /2/9^] where p = (n — r2 

{H)3g is 970 MeV in the Model I and 1311 MeV in the 



The wavefunction has the form 

/V2 and A = (n + fa - 2f3) /v^. 
The minimum for rriN 
Model II respectively, reached at the same f3 in both models. The same value of 
f3 is also used for the 6q system. This is equivalent with imposing the "stability 
condition" which is of crucial importance in resonating group method (RGM) cal- 
culations [1,2]. The quantity K^ei represents the relative kinetic energy of two 3^ 
clusters separated at infinity 



K, 



rel 






(20) 



where m above and in the following designates the mass of the u or d quark. For 
our value of P this gives Krei = 0.448 GeV. 



Cluster model 

At Z ^ oo the symmetries corresponding to baryon-baryon channels, namely 
[51)^5 and [SSji^s, must appear with proper coefficients, as given by Eq. (21). 
The contribution due to these symmetries must be identical to the contribution 
of V^ to two nucleon masses also calculated with the Hamiltonian (11). This is 
indeed the case. In the total Hamiltonian the contribution of the [411]j7'5 V^ state 
tends to infinity when Z ^ oo. Then this state decouples from the rest which is 
natural because it does not correspond to an asymptotic baryon-baryon channel. 
It plays a role at small Z but at large Z its amplitude in the A^A^ wavefunction 
vanishes, similarly to the "hidden colour" states. Actually, in diagonalizing the 
total Hamiltonian in the basis (2)- (5) we obtain an A^A^ wavefunction which in the 
limit Z — * oo becomes [14] 



^ 



NN 



1 

3 



momps) + o l[42]o[33]^5) - - |[42]o[51]^^) 



(21) 



The adiabatic potential drawn in Figs. 1 and 2 is defined according to Eq. (19) 
where {H)z is the lowest eigenvalue resulting from the diagonalization. Fig. 1 
corresponds to 5* = 1, / = in the Model I and Fig. 2 to the same channel in the 
Model II. Note that from these curves one should subtract Krei of Eq. (20) in order 
to obtain the asymptotic value zero for the potential. Similar results are obtained 
in the SI = (01) channel. One can see that the potential is repulsive at any Z in 
the Model I, but a small attractive pocket appears in the potential of the Model II 
(see the zoom in the inside box of Fig. 2). 




FIGURE 1. Comparison of the adiabatic potential in Model I, for SI = (10), calculated in the 
cluster model (solid curve) and the molecular orbital basis (dashed curve). 




FIGURE 2. Same as Fig. 1 but for the Model II. The upper box represents a zoom of the small 
attractive pocket. 

Molecular orbital basis 

In the molecular orbital basis the asymptotic form of ipNN is also given by Eq. 
(21) inasmuch as r ^ /? and I ^ L as indicated below Eq. (10). 

We diagonalize also the Hamiltonian and use its lowest eigenvalue to obtain the 
A^A^ potential according to the definition (19). The S* = 1, / = case is illustrated 
for Models I and II in Figs. 1 and 2 respectively, for a comparison with the cluster 
model basis. The repulsion reduces by about 15 % in the ^5*1 channels when passing 
from the cluster model basis to the molecular orbital basis. From Figs. 1 and 2 
one can see that the molecular orbital basis has an important effect up to about 
Z K, 1.5 fm giving a lower potential at small values of Z. For Z k, \ hn. it gives 
a potential larger by few tens of MeV than the cluster model potential. However 
there is no attraction at all in the Model I (for details see Ref. [18]). 



MIDDLE RANGE ATTRACTION 



In principle we expected some attraction at large Z due to the presence of the 
Yukawa potential tail in Eq. (14). This is indeed the case in Model II. In Model 
I, we adopt a consistent procedure assuming that besides the pseudoscalar meson 
exchange interaction there exists an additional scalar, a-meson exchange interaction 
between quarks. This is in the spirit of the spontaneous chiral symmetry breaking 
mechanism on which the GBE model is based. The a-meson is the chiral partner 
of the pion and it should be considered explicitly. 



Actually once the one-pion exchange interaction between quarks is admitted, one 
can inquire about the role of at least two-pion exchanges. Recently it was found 
[21] that the two-pion exchange also plays a significant role in the quark-quark 
interaction. It enhances the effect of the isospin dependent spin-spin component 
of the one-pion exchange interaction and cancels out its tensor component. Apart 
from that it gives rise to a spin independent central component, which averaged over 
the isospin wave function of the nucleon it produces an attractive spin independent 
interaction. These findings also support the introduction of a scalar (o"-meson) 
exchange interaction between quarks as an approximate description of the two- 
pion exchange loops. 

For consistency with the parametrization [8] we consider here a scalar quark- 
quark interaction of the form 



VJr) 



9l 



1 



A-K 12mimj 



{e{r - rd)^^ 



-fJ^a-r 
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(22) 



where n^ = 675 MeV and Tq, a' and the coupling constant 5'^/47r are arbitrary 
parameters. In order to be effective at medium-range separation between nucleons 
we expect this interaction to have Tq ^ tq and a' ^ a. Note that the factor l/niimj 
has only been introduced for dimensional reasons. 

We first looked at the baryon spectrum with the same variational parameters 
as before. The only modification is a shift of the whole spectrum which would 
correspond to taking Vq ~ —60 MeV in Eq. (12). 

For the 6g system we performed calculations in the molecular basis, which is more 
appropriate than the cluster model basis. We found that the resulting adiabatic 



0.6 



0.3 



Z(fm) 



FIGURE 3. The adiabatic potential Model I in the molecular orbital basis for SI = (10) (full 
curve) and SI = (01) (dashed curve) with pseudoscalar + scalar quark-quark interaction. 



potential is practically insensitive to changes in n„ and Tq but very sensitive to a'. 
In Fig. 3 we show results for 

r', = 0.86 fm, a' = 1.47 fm'\ gl/An = gl/Aix (23) 

One can see that V^ produces indeed an attractive pocket, deeper for SI = (10) 
than for (01), as it should be for the A^A^ problem. The depth of the attraction 
depends essentially on a'. The precise values of the parameters entering Eq. (22) 
should be determined in further RGM calculations. As mentioned above the Born- 
Oppenheimer potential is in fact the diagonal RGM kernel. It is interesting that 
an attractive pocket is seen in this kernel when a a-meson exchange interaction is 
combined with pseudoscalar meson exchange and OGE interactions (hybrid model), 
the whole being fitted to the A^A^ problem. The same interaction should also be 
introduced in Model II, but here, our purpose is to show how the a-exchange 
interaction gives some middle-range attraction. 

SUMMARY 

We have calculated the A^A^ potential in the adiabatic approximation as a func- 
tion of Z, the separation distance between the centres of the two 3g clusters. We 
used two different parametrizations of a constituent quark model where quarks in- 
teract via pseudoscalar meson exchange. The orbital part of the six-quark states 
was constructed either from cluster model or molecular orbital single particle states. 
The latter are more realistic, having the proper axially and reflect ionally symme- 
tries. Also technically they are more convenient. We explicitly showed that they 
are important at small values of Z. In particular we found that the A^A^ potential 
obtained in the molecular orbital basis has a less repulsive core than the one ob- 
tained in the cluster model basis. However none of the bases leads to an attractive 
pocket in one of the parametrizations considered here. We have simulated this 
attraction by introducing a a-meson exchange interaction between quarks. 

The present calculations give us an idea about the size and shape of the hard 
core produced by the GBE interaction. Except for small values of Z the two bases 
give rather similar potentials. Taking Z as a generator coordinate the following 
step is to perform a dynamical study based on the resonating group method which 
will provide phase-shifts to be compared to the experiment. These calculations are 
underway. 
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